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We study the impact of TeV-scale sterile neutrinos on electro-weak precision observables and 
lepton number and flavour violating decays in the framework of a type-I see-saw extension of the 
Standard Model. At tree level sterile neutrinos manifest themselves via non-unitarity of the PMNS 
matrix and at one-loop level they modify the oblique radiative corrections. We derive explicit 
formulae for the S, T, U parameters in terms of the neutrino masses and mixings and perform a 
numerical fit to the electro- weak observables. We find regions of parameter space with a sizable 
active-sterile mixing which provide a better over-all fit compared to the case where the mixing is 
negligible. Specifically we find improvements of the invisible Z-decay width, the charged-to-neutral- 
current ratio for neutrino scattering experiments and of the deviation of the W boson mass from 
the theoretical expectation. 
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I. INTRODUCTION 

The Standard Model (SM) is extremely successful and 
has passed numerous experimental tests. Moreover the 
last missing piece, the Higgs particle, has recently likely 
been seen by the ATLAS and CMS collaborations [1 . 
On the other hand, the SM is for a number of theoret- 
ical reasons incomplete. It also docs not explain three 
known experimental facts, namely the tiny active neu- 
trino masses, the baryon asymmetry of the Universe and 
the existence of dark matter. A simple yet elegant way to 
solve two or even all of these problems is to supplement 
the SM by three singlets: 



(1) 



Nf are Majorana fields, £q, are the lepton 
i(j2(j)* is the conjugate of the Higgs doublet. 



where N, 
doublets, 

and h are the corresponding Yukawa couplings. 

This is a small modification of the SM, also because so 
far the principles which select the fermionic representa- 
tions in a given theory are (besides for anomaly condi- 
tions) unknown. Nevertheless, this modification allows to 
accommodate those three experimental facts: First, after 
breaking of the electro- weak symmetry the active neutri- 
nos acquire masses via the typc-I sec-saw mechanism 0- 
1^. Second, within the same setup the observed baryon 
asymmetry of the Universe can be naturally explained 
by leptogcnesis Q • According to this scenario the lepton 
asymmetry which is produced by the heavy neutrinos 
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is converted into the bar yon asymmetry by anomalous 
electro- weak processes |8l-[l(]l|. Depending on the values 
of the masses and couplings in ([T]) the generation of the 
asymmetry proceeds either through decays [ll| - |23| or via 
oscillations of the heavy neutrinos [1^. The extension 
of the SM spectrum by sterile neutrino states may even 
provide a solution to the third problem, since sterile neu- 
trinos are a perfect dark matter candidate (see e.g. f25l - 
[3^). A further possible implication of sterile neutrinos 
is a modified active neutrino flux in reactor experiments 
[33| . This could explain the observed reduced electron 
antineutrino fluxes [34| . 

Note that the addition of singlets to the SM has actu- 
ally important consequences. Not only arc there many 
new parameters, but also the global symmetries are 
changed, since lepton number is broken. Furthermore the 
single scale of the SM, namely the electro-weak VEV, is 
amended by the new mass scales Mi. Within the type-I 
see-saw the sterile neutrinos pick up masses proportional 
to Mi and these could in principle have any value which 
is not excluded by experiment. It is therefore interest- 
ing to study implications of sterile neutrinos assuming 
any mass. Some effects exist even for ultra heavy sterile 
neutrinos, while some are only phcnomcnologically im- 
portant for not so heavy states, which are in the mass 
range of several TeV or lower. 

The existence of Majorana neutrinos has well known 
consequences on the phenomenology below the electro- 
weak scale. In particular, the new states can contribute 
to the amplitude of the neutrinoless double-beta d ecay 
[35l - [39| and induce rare charged lepton decays [10, l4l| . 
Furthermore, they can affect the electro-weak precision 
observables (EWPOs) via tree-level as well as loop con- 
tributions and thus provide an explanation for anomalies 
in the experimental data. In particular, the tree-level 
effects result in non-unitarity of the active neutrino mix- 
ing matrix j40l | and lead to a suppression of the invisi- 
ble Z-decay width. This is in agreement with the long 
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standing fact that the LEP measurement of the invisible 
Z-decay width is two sigma below the value expected in 
the SM [i^l ■ Furthermore the neutral-to-charged-current 
ratio in neutrino scattering experiments can be changed 
thus providing an explanation for the NuTeV anomaly 
[4^. Also a slight shift of the W boson mass from the 
value derived from other SM parameters is induced, re- 
ducing the tension between the input parameters of the 
electro-weak fit and the experimentally observed value 



II. 



OBSERVABLES 



44| . These electro- weak observables are affected not only 
at tree-level but also by loop effects. The latter can be 
taken into account in the form of the oblique corrections 
and can partially 'screen' the tree-level contributions, as 
we will see later. The impact of oblique corrections on 
models with non-universal neutrino gauge coupling has 
been studied in j45l - i47j . In these works the tree- level ef- 
fects (ce, and Cr) and oblique {S, T and U) corrections 
have been treated as independent parameters. However, 
in the model described by the Lagrangian ([IJ they are 
functions of the Majorana masses and Yukawa couplings 
and are therefore not independent. 



Encouraged by the fact that sterile neutrinos arc very 
well motivated we study their phenomcnological impact 
in this paper. Specifically we consider TcV-scalc sterile 
neutrinos with a sizable active-sterile mixing and deter- 
mine their over-all contributions to the EWPOs and to 
indirect detection experiments in the framework of the 
see-saw type-I extension of the SM. In contrast to the 
previous studies we derive expressions for the T, U- 
parameters in terms of the masses and mixing angles 
of the Majorana neutrinos, thus all quantities effecting 
the observables are functions of the parameters in ([l}. 
Our approach is therefore completely self-consistent. In 
Sec.|ll]we discuss the influence of heavy Majorana neutri- 
nos on the phenomenology below the electro-weak scale 
and provide expressions for the S*, T, [/-parameters. Next 
we perform in Scc. lIIIl a likelihood fit to the EWPOs tak- 
ing into account up-to-date experimental constraints in- 
cluding neutrino oscillation data and limits from the non- 
observation of rare charged lepton decays and neutrino- 
less double-beta decay. We find that TeV-scale sterile 
neutrinos improve the overall fit by bringing the invisi- 
ble Z-decay width, the charged-to-neutral current ratio 
for neutrino scattering and the W boson mass in agree- 
ment with the experimentally observed values within 
the experimental precision. The best-fit regions provide 
testable experimental signatures. For the normal- and 
quasi-degenerate light neutrino mass spectra we find that 
decay rates are close to the current experimental 
sensitivities. For the inverted hierarchy the mass region 
is such that part of the parameter space can be tested at 
the LHC after the 14 TcV upgrade. Finally, in Sec.HVl 
we summarize the obtained results and present our con- 
clusions. 



Despite the fact that the Majorana neutrinos in ([TJ are 
SM singlets, after the electro- weak symmetry breaking 
they mix with the active neutrinos through the induced 
Dirac mass terms and therefore also couple to the Z and 
W bosons. Expressed in terms of the mass eigenstates 
the corresponding part of the Lagrangian takes the form 



W,Y.Zli:c.^e,,,r^^^LrPLe^ + h.C. , (2) 



where e„ denote the charged leptons, Ui denote the light 
(for z < 3) as well as heavy (for 4 < z < 3 + n) neutrinos, 
and U is the full unitary (3 -I- n) x (3 -I- n) neutrino mixing 
matrix. To give masses to at least two active neutrinos, 
as required by oscillation experiments, we need two or 
more heavy neutrinos, i.e. n > 2, but otherwise n is un- 
constrained. Below we review the phenomenology of ^ 
together with the corresponding up-to-date experimental 
results. 

Lepton-flavor violating decays. In the scenario under 
consideration the branching ratio of ^ — > 67 decay is gi- 
ven by Hi] 

BR(^ ^ 67) = — •^-^—IS.W (3) 

i (/i — > eiyiy) SZir 

where 6„ = 2 J^i UgiU^; g (mf/M^) and the loop func- 
tion g is defined by 

= f Ji^^l^[2(l - a)(2 -a) + a{l + a)x]. 

Note that we use m,; to denote the masses of both light 
and heavy neutrinos. Since the masses of the active neu- 
trinos are very small we can neglect them in the loop 
integral, g{fnf /M^^) « g{0) ~ 5/3. Using unitarity of 
the full mixing matrix we then find 



= 2E-=;u:,;U^. [g (mj/M^) - 5/3] 



(4) 



The recent limit on this branching ratio obtained by the 
MEG collaboration El is 



BR(/_t^ 



^7) < 2.4-10" 



(5) 



at 90% confidence level. From Q and (O we can infer 
bounds on the products of the mixing elements U^^ and 
Uei. An analogous relation also exists for the t — > 67 
decay. However, the corresponding experimental con- 
straints are much weaker and will not be considered here. 

Neutrinoless double-beta decay. Neutrinoless double- 
beta decay constrains the effective mass of the electron 
neutrino (mee). The latter receives contributions from 
the light as well as from the heavy mass eigenstates [s^ : 



l('^ee)|«|EllU^ 



(6) 
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For masses of the heavy neutrinos in the TeV range one 
can use an approximation: F{A,mi) « {ma/mi)^ f[A), 
where ttIq k, 0.9 GeV and f{A) depends on the decay- 
ing isotope under consideration [39l [soj . A conservative 
bound, |(TOee)| < 0.4 eV, has been recently obtained by 
EXO collaboration [U. 

Unitarity and lepton universality violation. For the 
following analysis it is convenient to represent U in the 
form 



U 



(7) 



The [n x 3) matrix M describes the active-sterile mix- 
ing. An obvious consequence of a nonzero active-sterile 
mixing is that the (3 x 3) PMNS matrix is no longer 
exactly unitary (40| . The deviation from unitarity can be 
parameterized by 



Ei>4l''--'^ai| 



(8) 



In general the quantities eg, and t-r are not equal. In 
other words, there is also a violation of lepton universal- 
ity. Experimental bounds on linear combinations of the 
Eq read [43 



ee -e^ = 0.0022 ± 0.0025, 
e^-€r =0.0017 ±0.0038, 
e^-e^ = 0.0039 ±0.0040. 



(9a) 
(9b) 
(9c) 



The stringent experimental bound on the — > 67 branch- 
ing ratio implies that either or is negligible in 

Heavy neutrinos at colliders. Heavy singlet Majorana 
fcrmions have been searched for at the LEP and LHC 
colliders, see [s^, H^l for a review. The searches were 
based on the production of the heavy state due to a con- 
siderable mixing to the active neutrinos. 

At LEP the heavy neutrino could be produced in e+e^ 
annihilation, e~^e~ — ^ Nv, via s-channel Z-exchange as 
well as via t-channel W^-exchange. The produced Ma- 
jorana neutrinos then rapidly decay via the weak neu- 
tral or charged currents: N — >■ Zv and N — )• We. A 
search for heavy neutrinos with masses up to 200 GeV 
has been performed by the L3 and Delphi collaborations 
using the latter decay channel with W decaying into 
hadrons [13, HI]- The experimental signature of these 
events would be one isolated electron plus hadronic jets. 

New limits on the active-sterile mixing for the heavy 
neutrino masses up to 210 GeV have been obtained re- 
cently by the CMS collaboration [E^ using a dilepton de- 
cay channel with two leptons of equal charge and flavour 
plus jets, see Fig.[T] Violation of lepton number in this 
process occurs due to the Majorana nature of the sterile 
neutrino. For large Majorana masses the square of the 
momentum transfer in the propagator of the intermedi- 
ate neutrino can be neglected and the production cross 
section depends on the combination | U^j m~^\. As- 
suming that only one of the heavy neutrinos couples to 
the electron with a strength |Ueip « 5.2 -10"^ it has been 




FIG. 1: Lepton number violating process mediated by the 
sterile neutrino. 



estimated in [531 that after the LHC upgrade to -^s ~ 14 
TeV this channel can be used to search for Majorana 
neutrinos with masses up to roughly 800 GeV. For these 



parameters we get | J2i 



6.5-10-3 TeV^i and 



one event per 100 fb~^ is expected. Even though this is 
a very weak statistical signal it has been noted that in 
the highest mass region (800 GeV) the leptons are emit- 
ted back to back providing a very clean signature and 
allowing for an excellent background suppression. If the 
luminocity of 400 fb~^ is reached one can hope to find 
the heavy neutrinos even if the above defined combina- 
tion of the mixings and masses is as small as 3.25 • 10^'^ 
TeV-i. 

A possible way of testing even higher mass ranges is 
to study electron-proton collisions where the dominant 
detection channel is e~ + q — >■ q' + e'^ + W with the 
Majorana neutrino as an intermediate state. In (5^ it 
has been estimated that at an electron-positron collider 
with a center of mass energy sa 6 TeV mass ranges 
up to 1.5 TeV can be tested. 

Non-unitarity in neutrino oscillations. The standard 
oscillation formula assumes unitarity of the PMNS ma- 
trix which is violated in the model under considera- 
tion. The non-unitarity of ^ will manifest in a modified 
oscillation probability formula (40| . which in particular 
includes a zero- length effect [s^ : 



P^p{L = 0) = 



5al3 (1 - 2£a) + tgefj 

{l-eg){l-ep) 



(10) 



Being of second order in e for a ^ j3 this effect is strongly 
suppressed. As discussed by a number of authors, see [40[ 
and references therein, the probabilities do not add up 
to unity. At non-zero distances the non-unitarity effects 
are linear in e [60| and may affect the results of global 
fits. For the purposes of this work we will use the best- 
fit results for the light neutrino parameters obtained in 
a unitary fit as an input, see [6l[. The order e correc- 
tions to the parameters of the active neutrinos will be 
computed using (|20p . At the end, when comparing the 
corrected values to the results of the unitary fit in [6l| 
we find agreement within the one sigma intervals. The 
values are also consistent with the one sigma intervals 
obtained in a non-unitary fit [40| . Furthermore, we have 
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checked that our results are not sensitive to changes of 
the active neutrino parameters within their experimental 
errors, which justifies our procedure. 

Electroweak precision observables. An important con- 
sequence of non-unitarity is that the couplings of the light 
neutrinos to the Z and W bosons are suppressed with re- 
spect to their SM values. Taking this effect into account 
we find that the invisible Z-decay width is suppressed: 



Tinv/ [rinvlsM ~ 3 '^ai^ ^a) 



(11) 



This means that the 'effective number of neutrinos' 
measured in the .Z-decay is slightly less than three, 
which is qualitatively agreeing with the LEP results 
[4^ . Similarly, cross-sections of the charged- and neutral- 



EWPO 


Theory (Standard Model) 


Experiment 


Tiept (MeV) 


84.005 ± 0.015 


83.984 ± 0.086 


Tinv/Tiept 


5.9721 ± 0.0002 


5.942 ±0.016 


sin^ 


0.23150 ± 0.0001 


0.2324 ± 0.0012 


gl 


0.3040 ± 0.0002 


0.3026 ± 0.0012 




0.0300 ± 0.0002 


0.0303 ± 0.0010 


Mw (GeV) 


80.359 ±0.011 


80.385 ± 0.015 



TABLE I: Theoretical predictions and experimental results 
for electro- weak precision observables (EWPO). The theoret- 
ical predictions are taken from [g^l. The experimental values 
for the invisible and leptonic Z-dec&y widths as well as the 
Weinberg angle are from For the Weinberg angle we use 
the value measured in the hadronic processes to make sure 
that it is free of the non-unitarity corrections. The values of 
ql and qr are taken from the NuTeV results after including 
a recent NNLO analysis [i^ . [63 | . 

current neutrino scatterings on quarks are also affected: 



cc 



NC ^ \2 



(12a) 
(12b) 



The stronger relative suppression of the neutral current 
interactions can be observed in experiments measuring 
ratios of the corresponding cross sections. In particu- 
lar, it is qualitativelyconsistent with the results of the 
NuTcV experiment [6J-[63]. Another important conse- 
quence of non-unitarity is that - the Fermi con- 
stant measured in the muon decay - is not equal to 
the Fermi constant measured in experiments with semi- 
leptonic processes, but 



(13) 



Since the muon decay width is used as an input in the SM 
fits, this modification influences many observables and 
has been used in to obtain bounds on the active- 
sterile mixing. However, as argued in j45l l46j . the im- 
pact of the heavy neutrinos is not limited to the above 
discussed tree-level effects. Heavy neutrinos contribute 
to the self-energies of the W and Z bosons and therefore 



modify their propagators. These loop effects can be de- 
scribed in terms of the S, T, U parameters [g^ . A more 
detailed description of the oblique correction formalism 
is presented in Appendix|Xl while the calculations are 
shown in Appendix|B] Combining the tree-level and one- 
loop contributions one obtains [45|, [ioj 



Icpt 



[ricptlsM 

Tinv/ricpt 

cptlsM 



[Finv/Fl, 



sm' 



2 ^Icpt 



[sin^ e':^'] 



SM 



9l 



SM 



9l 



[SflJsM 



w 



WJsM 



: 1 ± 0.6 (ee±eM + 0-0145 T) 
0.00215", 

: 1 -0.67(ee±e^±er) 

0.00215"- 0.0015 r, 
: 1-0.72 (ee±e^± 0.0145 T) 

0.0016 5", 

: 1 ± 0.41 Ee- 0.59 

0.0090 5 ± 0.0022 r, 
:l-1.4ee-2.4e^ 

0.0315-0.0067T, 
: 1 ± 0.11 ee± 0.11 

0.0036 5 ± 0.0056 r± 0.0042 C/. (14f) 



(14a) 



(14b) 



(14c) 



(14d) 



(14e) 



Comparing pTjl and (|14b[) we see that the first line in 
the latter equation is the leading-order expansion of the 
former one and that the loop corrections enter through 
the T and 5" parameters. 

Cancellation mechanism. Since 5" and U are related to 
derivatives of the W and Z boson self-energies whereas T 
is proportional to a difference of the two, usually 5" and 
U arc subdominant in comparison to T. If 5 is neglected 
then (|14ap and (|14cp depend on the same combination 
of tree-level and one-loop corrections, ee ± e^^ ± 2aemT, 
where aem is the fine structure constant. The reason is 
that in this approximation the shift of these observables 
is solely due to a shift in [i^ . It has been argued in 
[40| that a shift in G^ has dramatic effect on the electro- 
weak observables and thus Ce and are strongly con- 
strained. However, as follows from ([H]). if the tree- level 
and radiative contributions are of a similar size the shift 
induced by the tree-level effects can be 'screened' by a 
sizable negative T parameter. Typically, most of the SM 
extensions result in a positive T parameter. Therefore, it 
was assumed in [46| that the Higgs mass is much larger 
than the used reference value, niH = 115 GeV. However, 
interpreting the recent discovery of the ATLAS and CMS 
collaborations as the SM Higgs, this option is excluded. 

Majorana neutrino contribution to the S, T, U param- 
eters. Here we argue that the screening can be realized 
without invoking any new physics beyond that already 
introduced in ([ij. As demonstrated in [69| for a model 
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with a single heavy neutrino, in a certain range of the 
parameter space Majorana neutrinos can generate a neg- 
ative contribution to the T parameter. In the model with 
n heavy neutrinos we find 



parameter we obtain 

Tt 



Tm + T, 



1 



SM 



+ Eo^^lMO, ml,ml) - EMMl ml,ml) 
- 2(M2/A/w)'EaULu„,AQ(A'/2^,m2,m2)] , 



X [E-t:iE.,3ULu„,.Ut^U^,Q(0,m?,m^2) 

+ E ■,l:iEa^ uLu„,-Ut^U,,m,m,-i?o(0, mlmj) 

-2E-=i"E.uLUo.g(0,m2,m2) 

+ Ea^lBo{0,<,0], (15) 

where to^ denote masses of the charged leptons. To shor- 
ten the notation we have introduced 

Qiq^,ml ml) = {D - 2)B22{q\ml ml) 

+ q^[Bi{q\mlml) + B2i{q\mlml)] , (16) 

where B^, Bi, B21 and B22 arc the usual one- loop func- 
tions [zOl, D = 4 — 2e and e — > 0. Details of the calcula- 
tion can be found in Appcndix|B] Note that although 
the loop functions are divergent and contain an arbi- 
trary scale /i, their combination (jlSp is finite and inde- 
pendent of ^. This can be shown using the unitarity of 
the full mixing matrix and the type-I see-saw condition 



,m2) 



(19) 



(mi) 



For a single generation Eq s. (1151) . (flS)) and (fT9|) reduce to 
the expressions derived in [69| . The observation that the 
heavy neutrinos can 'screen' the tree-level contributions 
by generating a negative T parameter together with the 
explicit formulae for the S, T, U constitutes one of the 
main results of the present work. 



III. FIT TO THE OBSERVABLES 

In this work we investigate the impact of TeV-scale 
sterile neutrinos with a sizable active-sterile mixing on 
the physics below the electro- weak scale in the framework 
of the see-saw type-I extension of the SM. In the see-saw 
limit the active-sterile mixing is of order rh^j {Mfi)~^ . On 
the other hand, the see-saw formula for the mass ma- 



Ei=r Ua»"iiU,,^ = with a,/3 e {e,//,r}. trix of the light neutrinos reads = -mJ){MR)-^mD- 



Therefore, to compute (jlSp one can simply drop the 
terms in the expansion of the loop functions and eval- 
uate the remaining integrals numerically or analytically. 
As indicated in (|T5l) . the T parameter is a sum of the 
SM and new contributions. The SM contribution can be 
obtained from psp by setting the active neutrino masses 
to zero and taking into account that = and ^ is 
unitary in the SM: 



TsM — — 



1 



[3O(0,0,0)-2E„Q(0,0,m2) 
+ E^<Bo{0: <:<)]■ (17) 



This expression is also finite. Note that the PMNS ma- 
trix docs not appear in (|17p since the mixing becomes un- 
physical for massless neutrinos. The S parameter reads 



Comparing the two expressions we conclude that for a 
large active-sterile mixing the smallness of the active neu- 
trino masses can not be explained by the scale suppres- 
sion. Instead, contributions from different heavy states 
have to mutually cancel. To ensure the cancellation more 
than one sterile neutrino is needed. As already men- 
tioned above, two heavy neutrinos are required to ensure 
that at least two of the light neutrinos are massive. We 
have also seen that the mixing to either the electron or 
the muon neutrino has to be strongly suppressed to be 
compatible with the current ^ — >■ 67 bounds. Combined 
together, the two conditions imply that the cancellation 
is only possible for 71 > 3. We consider n = 3 in this 
paper, which arises naturally in the framework of many 
left-right symmetric or GUT models. 

For n = 3 both ^ and ^ are (3 x 3) matrices which 
can be written in the form iTlll: 



1 



Stot ^ + SsM = - ^^^^2 

+ Ea"^' ^0(0, ml,ml) + EaQiMl mlml) 
-2j:^mlBo{Mlml,ml)], (18) 

where /\Q{q^ ,m\,ml) = Q{0,m\,m\) — Q{q^ ,m\,mD 
and ASq is defined in the same way. The SM contribu- 
tion, SsM, is calculated analogously to Tsm- For the U 



heavy : 



^ = {I - 3£ 3?.^Y U 



(20a) 
(20b) 



with O being an arbitrary complex orthogonal matrix, lA 
the unitary matrix diagonalizing m^, Whoavy the diagonal 
mass matrix of the heavy neutrinos and mught that of 
the light neutrinos. The entries of lA are extracted from 
a global fit to atmospheric, reactor and solar neutrino 
oscillation data. Note that a possible non-unitarity of 
is neglected in usual global fits, i.e. it is assumed that 
= U. The resulting best fit values and la ranges for 
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the three mixing angles and the Dirac CP-phase^ read 
(see e.g. [6l|) 



where BR*''(/Lt -J> 67) and 



are the theoretical pre- 



sin^ 6*12 
sin^ 6*23 
sin^ 6*13 



JCP 



0.30 ±0.013, 

n 41+0.037 
^■^^-0.025 ' 

0.023 ±0.0023, 



300 



+66 
-138 ' 



(21a) 
(21b) 
(21c) 
(21d) 



whereas the two Majorana phases arc unknown and will 
be set to zero in the fit. The remaining degrees of freedom 
are the three masses of the heavy neutrinos and the en- 
tries of the matrix O. The latter can be parametrized by 
three complex angles. The freedom of choosing O may be 
used to suppress the active-sterile mixing for one active 
neutrino flavour. Whether this suppression is possible or 
not depends on the structure of the light neutrino mass 
matrix, as can be inferred from (1201). Since the absolute 





NH 


IH 


QD 


mi (eV) 


~ 


4.85 ■ 10"^ 


~ 0.1 


m2 (eV) 


8.660 ■ 10-=^ 


4.93 ■ 10"^ 


~ 0.1 


ma (eV) 


4.97- 10"^ 


~ 


~ 0.1 



TABLE II: Assumed masses of the light neutrinos for the nor- 
mal (NH), inverted (IH) and quasi-degenerate (QD) neutrino 
mass spectra. 

neutrino mass scale is unknown at present, we consider 
three possible neutrino mass spectra called NH, IH and 
QD which are defined in Tablellll In the case of the QD 
mass spectrum, the masses are constrained by the recent 
WMAP bound < 0.44 eV [zl. For the numerical 

example in the QD case we use rrii ~ 0.1 eV. 

The minimization of the function is performed with 
a statistical method described in Appendix|Dj The gen- 
eral approach is as follows: For a given point in the pa- 
rameter space we compute the — >■ 67 branching ratio 
and the effective electron neutrino mass {niee), relevant 
for the Qvpp. For most of the points with a sizable active- 
sterile mixing one of these bounds is violated. To take 
the two bounds into account in the analysis we define 



2 ^ / BR*^(/i^e7)-BR-^(/.^e7) \ 
^"^"^ ~ y BR'="P(^ ^ 67) j 

X 6'(BR*''(^ 67) - BR'=^P(/.i 67)) , (22a) 
2 _ i\{K'l)\-\{ml7)V^ 



x0i\{ml':)\^\{mT/)\)., 



^ The one sigma best fit region for the CP-phase was determined 
in [T^ I to be 0.77 tt — 1.36 tt. However, the choice of the CP-phase 
has no significant impact on our results. 



dictions computed at the chosen point of the parameter 
space, whereas BR*^'''^(/i 67) and |(mggP)| are the cor- 
responding Icr experimental upper bounds. We use the 
theta step function to restrict the contributions to the 
total to cases when the theoretical prediction exceeds 
the one sigma exclusion limit. Additionally, we check 
whether the universality constraints, see Eq. are ful- 
filled. Finally, we compute the S, T, U parameters and 
the corrected values of the electro- weak precision observ- 
ables Oi , see Eq. dH]) . The corresponding Xewpo value 
is calculated using 



Xewpo 



E 



(23) 



where Oi^sM denotes the predictions of the SM, 50i_sM 
are the theoretical errors and 50i are the experimental 
errors, see TablelH Note that we neglect off-diagonal el- 
ements of the covariance matrix in these contributions. 
The total x^ is given by the sum of ^ and ((23)) . 

Using this definition we find that for points with negli- 
gibly small active-sterile mixing ('natural' see-saw) ~ 
Xo = 7.5 (7.6 for the QD mass spectrum). This relatively 
large value of x^ is induced primarily by the anomalies 
in the invisible Z-decay width, the NuTeV results and 
the deviation of the W boson mass from the SM expecta- 
tion. The slightly higher value for the QD mass spectrum 
stems from the fact that neutrino masses of the order of 
0.1 eV induce an effective electron neutrino mass (ml^^) 
comparable to the current upper bound. 

To get a rough estimate of the goodness of the fit we 
compute the ratio of x^ to the number of degrees of free- 
dom. The S and U parameters are always very small 
and can be neglected in the fit. Therefore, the initial set 
of free parameters, the three masses and three complex 
angles, maps to four quantities: Ce, e^, er and T. Since 
we fit six observables and two constraints, the number of 
degrees of freedom is (6 + 2) — 4 = 4 and Xo/d.o.f « 1.9 
for the case of negligible active-sterile mixing. Below we 
perform a fit for each of the three light neutrino mass 
spectra to all EWPOs and constraints. In the following 
in all plots the full x^ is presented. 

For all mass spectra we will show how the EWPOs in 
the best-fit regions are shifted with respect to the SM 
predictions. The question will be addressed whether a 
direct detection or an indirect signal of the new states 
is feasible in near future. To study this possibility the 
/i ^ 67 branching ratio, the 0/3^1/ decay rate and the 
strength of a collider signal will be estimated. 

Normal mass hierarchy. Let us first consider the nor- 
mal mass hierarchy of the light neutrinos. The /i — > 67 
(22b) branching ratio is suppressed in particular if is small. 

As discussed above, if S and U are neglected then the 
Z boson leptonic decay width (|14ap and the Weinberg 
angle (|14c|) depend on the same combination of tree-level 



and one-loop corrections, + + 2Q!e„jT. The reason 
is that in this approximation the shift of these observ- 
ables is solely due to a shift of with respect to Gi?-. 
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In Fig. [2] we show the fit to Fiopt and the Weinberg angle 
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FIG. 2: s-s a function of ee + ep. and 2aemTN (NH), here 
only Fiopt and the Weinberg angle are fitted. 

only. The plot shows the (approximately) linear relation 
between eg + and the T parameter. The fact that the 
band does not start at the origin indicates that the two 
observables favour a slightly negative T in general even 
in a different new physics model. This plot demonstrates 
how a negative value of the T parameter can screen the 
contributions of the neutrino mixing + e^x to the above 
mentioned observables. We will see that this cancella- 
tion mechanism is still at work when more observables 
are considered. Large cancellations are present when the 
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FIG. 3: for three d.o.f. as a function of te + e^ and 2aemljv 
(NH). Here the W boson mass is excluded from the fit. 



other EWPOs except the W boson mass are taken into 
account. However, in this case the simple linear depen- 
dence is violated due to the necessity of fitting the other 
observables, see Fig.[31 The W boson mass has been mea- 
sured with a very high accuracy and a large negative T 
would induce unacceptably large corrections to it. The 
fit excluding Mw is of interest since an extended scalar 
sector appearing, for instance, in the type-I-fll sce-saw 
extension of the SM will shift the W mass and might ab- 
sorb the additional negative T contribution restricted in 
the type-I see-saw. 
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FIG. 4: for four d.o.f. as a function of te -f and 2ae 
(NH), fit to all EWPOs in (fTI)). 



In Fig. m we show how the cancellation pattern chan- 
ges if all observables in ([T4)) are included. In this case for 
the best-fit points we have ~ 4.0 and a small negative 
T is favored. The latter partially compensates the small 
positive contribution of -t- to G^. 

As can be inferred from Fig.[6l at the best-fit points 
at least one the Majorana neutrinos is relatively light, 
around one TeV, and has a sizable coupling to at least 
one of the charged leptons, see Fig. [5] In other words, the 
current data favour the low-scale type-I see-saw with a 
considerable active-sterile mixing over the standard sce- 
nario ('natural' see-saw) where this mixing is negligible 
and the masses of the heavy neutrinos are close to the 
GUT scale. On the other hand, as can be seen in a 
different projection of the same plot in Fig. [SI heavy neu- 
trinos with masses below TeV scale are also disfavoured. 
The reason is that for sizable active-sterile mixing in this 
mass range experimental bounds are more stringent. As 
already mentioned above, for a sizable active-sterile mix- 
ing the cancellation of the contributions to the light neu- 
trino masses imposes constraints on the mass spectrum 
of the heavy ones. Due to requirements of suppressed 
Ee/^ an active-sterile mixing pattern occurs where the 
first and third heavy mass cigenstatcs have comparable. 
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FIG. 5: The lightest heavy neutrino mass as a function of 
for four d.o.f. and + Sfj, + tr (NH). Here is suppressed. 
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FIG. 6: Ee + e/j + ^T as a function of the lightest heavy neutrino 
mass and for four d.o.f. (NH). Here is suppressed. 



sizable mixing to the active flavours vifhile the mixing 
of the second is small. We find that the first and third 
heavy mass eigenstates have approximately equal masses 
whereas the mass of the remaining sterile neutrino is con- 
siderably larger. This mass pattern leads to a small neg- 
ative T parameter as discussed above. 

At the best-fit points the deviation of the PMNS ma- 
trix from unitarity is of the order of 



/ 1 - IQ- 



10-" 2 • 10-3 
1 • IQ-"* 1 - 10-5 2 - 10-4 
y 2- 10-3 2-10-'^ 



3 1 

-4 



which is in excellent agreement with the bounds derived 
in [t^ ]. This consistency check justifies our procedure 
and the applicability of the approximation in (|20p . The 
corresponding (largest) zero-length transition probability 
is P^^-^i,^ K, 7 • IQ-^. This value is too small to explain 
the current short-baseline anomalies (75l | but might con- 
tribute to the observed discrepancies. 

In Fig. [7] we demonstrate how the EWPOs in the con- 
sidered scenario are shifted from the values expected in 
the SM towards the experimentally observed ones. At 
the chosen best-fit point ~ 4.0 and therefore the ab- 
solute improvement of the fit is Ax^ ~ 3.5 leading to 
a x^/d-o.f « 1.0. Some of the observables are shifted 
away from the experimentally measured value compared 
to the SM prediction. Nevertheless, they remain within 
the one sigma interval around the experimental results. 
In addition, observables for which the SM prediction is 
in tension with data are brought into the one sigma in- 
tervals around their experimental values. This leads to a 
global picture in which all observables agree on the one 
(or in the case of g\ on the 1.2) sigma level with the 
experiments. This leads to a x^ /d.o.f of order one. The 
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FIG. 7: EWPOs calculated at the best-fit point for NH and 
suppressed e,j (green dots) compared to the experimentally 
observed values, denoted by the zero line. The coloured lines 
stand for the respective experimental sigma deviations, thus 
the displacement of the predicted values form the observations 
is presented in units of the experimental error. Note that for 
the Oi/fiP and /i — >■ 67 constraints we present only the one 
sigma exclusion limits. The theoretical predictions of the SM 
with their theoretical uncertainties, see TablelH are displayed 
as well (blue bars). (The best-fit point is at Mi — 20.3 TeV, 
M2 = 14.1 TeV, Ms = 21.0 TeV, = 2.1-10-^ = 3.0-10-^ 
and er = 4.5 • IQ-^) 

improvement stems primarily from er which leaves 
unaffected (and thus equal to Gp) but suppresses the 
invisible Z-dccay width. Furthermore, Ce shifts the W 
boson mass towards the measured value. The NuTeV 
observables are not significantly changed. The last two 
rows in Fig. [7] present the ratios of the induced effective 
electron neutrino mass {rriee) and the fj, ^ branch- 
ing ratio to the corresponding experimental bounds. It 
is interesting to note that for the normal mass hierarchy 
there arc many points in the parameters space with val- 
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FIG. 8: for four d.o.f. as a function of the ratios of the 
/I — >■ e7 branching ratio and |(mee)| to the corresponding 
experimental bounds (NH). Here e^j is suppressed. 



ues of these two quantities close to the threshold of the 
current experimental sensitivity, see Fig.|8l Knowing the 
masses of the heavy neutrinos and their couplings to the 
charged leptons we can also estimate the cross-section of 
the direct detection process depicted in Fig.[T] For all 
best-fit points the resulting value of | U^j "^f^l turns 
out to be at most of the order of 10~^ and thus beyond 
the LHC reach. 

The — >■ 67 branching ratio is also suppressed if eg 
is small. In this case we obtain a slightly worse best-fit 
value of « 4.8 which corresponds to x^/d.o.f w 1.2. 
As compared to the previous case, in the best-fit points 
the invisible Z-decay width is less suppressed and the 
shift in the W boson mass is smaller, sec Fig.|9l On the 
other hand, g\ is shifted towards the value measured by 
the NuTeV collaboration and thus brought in agreement 
with the data within the one sigma interval. 

Inverted mass hierarchy. For the inverted mass hier- 
archy scenario only the option of negligible can be 
realized and results in best fit values of ~ 5.5 and 
X^/d.o.f « 1.4. The EWPOs of the best-fit points are in- 
fluenced such that the predicted invisible Z-decay width 
is brought within the one sigma interval to the experi- 
mental value, see Fig.[Tni The W boson mass is shifted in 
the right direction but less significantly than in the nor- 
mal hierarchy scenario. On the other hand, is shifted 
towards the value measured by the NuTeV collaboration 
and brought in good agreement with the data. As can 
be inferred from Fig.[Tl] both the effective electron neu- 
trino mass (rriee) and the /i — >■ 67 branching ratio are in 
large parts of the parameter space close to the current 
experimental sensitivity. Furthermore, for a small frac- 
tion of the best-fit points, see Fig.[T21 there is a chance 
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FIG. 9: This is the analogue of Fig. [7] but for the case of 
suppressed ee in NH. At the best-fit point Mi = 4.1 TeV, 
M2 = 161.0 TeV, Ms = 7.1 TeV, = 1.0•10-^ = I.S-IQ-^ 
and er = 1.2 ■ 10"^. 



-3-2-10123 


Flept 

Finv/Flept 

sin^ O^f 
9l 
9% 
Ml 




• 




« 


1 



















FIG. 10: This is the analogue of Fig. [7] but for the case of 
IH and suppressed eg. At the best-fit point Mi — 551 GeV, 
M2 = 242 GeV, M3 = 469 GeV, te = 3•10"^ = 1.58- 10"^ 
and = 1.1 ■ 10"^. 



of observing the heavy Majorana neutrinos at the LHC 
after the planned upgrade to y/s = 14 TeV. 

Quasi- degenerate mass spectrum. For the quasi- 
degenerate mass spectrum it is quite difficult to find 
points with sizable that pass the — > 57 constraint. 
On the other hand, there are points with small and 
sizable Ce and Ct which satisfy the fi ^ e^, neutrino- 
less double-beta decay and universality constraints. For 
the best-fit points « 5 leading to an improvement of 
Ax^ « 2.6 and to x^/d.o.f « 1.25. At the best-fit points 
the shift of the invisible Z-decay width is smaller than 
in the normal hierarchy scenario, but still the prediction 
moves into the one sigma interval of the observations. 
The shift in the W boson brings prediction and exper- 
iment in excellent agreement, see Fig. 1131 The neutrino 
scattering observable is not significantly infiuenced in 
comparison to the standard model expectation. As can 
be inferred from Fig. [14] both the effective mass of the 
electron neutrino and — >■ 67 branching ratio here can be 
close to the exclusion limit. Note that the active neutri- 
nos induce a sizable contribution to the effective electron 
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FIG. 11: for four d.o.f. as a function of the ratios of 
the fi ^ ey branching ratio and |(mee)| to the corresponding 
experimental bounds in the case of IH and suppressed te. 
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FIG. 13: This is the analogue of Fig. [7] but for the case of 
the QD mass spectrum and suppressed e^. At the best-fit 
point Ml = 14.5 TeV, M2 = 10.5 TeV, M3 = 14.9 TeV, 
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FIG. 14: for four d.o.f. as a function of the ratios of the 
/X — >■ 67 branching ratio and |(mee)| to the corresponding ex- 
perimental bounds in case of quasi-degenerate mass spectrum 
and suppressed e^. 



FIG. 12: x^ for four d.o.f. as a function of mass of the heavy 
neutrino giving the leading contribution and coupling to the 
muon in the case of IH and suppressed te. The region above 
the solid line can be tested by the LHC after the upgrade to 
14 TeV if the anticipated luminosity ~ 400fb"^ is reached. 

neutrino mass |(?Tiee)| in this case and can constructively 
or destructively interfere with the heavy neutrino contri- 
butions depending on the Majorana phases of the active 
neutrinos. 

For all the best-fit points the resulting value of 
|^jU^jTO~^| turns out to be at most of the order of 
10~^ and thus beyond the LHC reach. 

It has been found that for all three mass spectra the 



current data favour a type-I see-saw with considerable 
active-sterile mixing over the a scenario where the mixing 
is negligible. For all three cases there are many points in 
the parameter space which generate signals in rare decay 
experiments close to the current sensitivities. For the IH 
a part of the parameter space can be tested at the LHC 
after the next upgrade. 



IV. SUMMARY AND OUTLOOK 

In this work we have studied the impact of TcV-scale 
sterile neutrinos with a sizable active-sterile mixing on 
electro-weak precision obscrvables and lepton number 
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and flavour violating decays in the framework of the type- 
I see-saw extension of the SM. The active-sterile mixing 
can in particular improve the global fit and reduce some 
anomalies in the experimental data: first, it results in 
a decreased value of the invisible Z boson decay width, 
which is preferred by the current data; second, it slightly 
increases the mass of the W boson and brings it in the 
one sigma interval around the observed value. Further- 
more, in the case when is not suppressed, the model 
affects the charged-to-neutral current ratio for neutrino 
scattering and brings gi^ in agreement with the values 
observed by the NuTeV experiment. 

In our analysis we have taken into account tree-level 
contributions of sterile neutrinos as well as their one-loop 
effects on gauge boson propagators. The tree-level con- 
tributions enter via non-unitarity of the PMNS mixing 
matrix, whereas the loop-contributions modify the prop- 
agators of the gauge bosons and can be taken into ac- 
count in form of the S*, T, JJ parameters. The oblique 
corrections can play an important role. In particular 
we have demonstrated that a sizable negative T param- 
eter 'screens' the effect of tree-level non-unitarity on the 
Fermi constant. This cancellation mechanism can recon- 
cile large active-sterile mixing with the current observa- 
tions. As compared to previous analyses we have derived 
explicit formulae for the oblique corrections in terms of 
the Majorana masses and mixings and therefore studied 
a full and self-consistent model. We have performed a nu- 
merical fit to the electro- weak precision observables tak- 
ing into account constraints from the non-observation of 
e^i and neutrinoless double-beta decay processes as 
well as constraints on lepton non-universality. Since the 
active neutrino mass hierarchy is unknown at present, we 
have considered the cases of normal, inverted and quasi- 
degenerate mass spectra. In all three cases regions of the 
parameter space with a sizable active-sterile mixing pro- 
vide a better overall fit to the data than regions where 
it is negligible. In other words, the current data favour 
the low-energy type-I see-saw with a considerable active- 
sterile mixing over the standard scenario ('natural' sec- 
saw) where this mixing is negligible. Together with the 
derivation of expressions for the oblique corrections in 
the type-I see-saw this finding is one of the main results 
of the present work. 

The is lowest for the normal hierarchy scenario and 
slightly higher for the inverted hierarchy and the case 
of the quasi-degenerate mass spectrum. We have also 
studied the experimental signatures of the best-fit points 
for each mass hierarchy. For all mass spectra the effec- 
tive electron neutrino mass (toec) ~ the quantity rele- 
vant for the neutrinoless double-beta decay - reaches at 
many best-fit points values that are close to the sensi- 
tivity threshold of current experiments. The expected 
active-sterile mixing is of the order of U^^y^^ ■ ~ lO^'^ in 
all three scenarios. The lightest sterile neutrino mass for 
normal- and quasi-degenerate mass spectra is around 2 
TeV. This ranges might become accessible at electron- 
proton colliders with energies above 6 TeV. For the in- 



verted hierarchy the Majorana masses can be as low as 
300 GeV. Thus, given the predicted mixing, this scenario 
can be partially tested at the LHC after the 14 TeV up- 
grade. 

The fit might improve in a type I-f II see-saw extension 
of the SM. since the VEV of the additional Higgs triplet 
also contributes to the W boson mass. Thus the negative 
T parameter can be larger and the cancellation mecha- 
nism more efficient. This extension will be studied in a 
forthcoming paper. 
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Appendix A: Oblique corrections 

The formalism of oblique corrections was developed 
to study new physics which affects the observables only 
via corrections to the propagators of the gauge bosons 
[tgI . Ft?! . The Lorentz structure of the corresponding self- 
energies is 

n^a?') = n,b(g2)gM- + {^f q^ terms) , (Al) 

where a, 6 = 7, , Z. Since in the considered processes 
all external fcrmions are light, the contributions of the 
q^(f terms is proportional to mj /My/,z and can safely be 
neglected. For the same reason, contributions of the box 
diagrams as well as the vertex corrections stemming from 
the new physics are negligible. Therefore, it is sufficient 
to consider only Tiab- The latter can be represented in 
the form 

n,,(g2)^nff (g2) + OT,,(g2), (A2) 

where {ab) = (77), (Z7), {ZZ),{WW). It follows from 
the Ward identity that (5n^^(0) = i5n^z(0) = which 
further constrains the number of independent quantities. 
It can be shown that for heavy new physics there are only 
three independent combinations of the above quantities 
that enter the expressions for the electro-weak observ- 
ables [7^. These are denoted by 5, T and U and read 
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Ml 



nzz(o) + n^^(Af|) 



Note that only for the T parameter the expression in 
terms of the unrenormahzed self-energies has the same 
form as in terms of the renormalized ones, compare ()A3bp 
and (|A7bp . 



nzz(o) iiww{Q) 



-nz^(M|) 



Ml 



1 UwwiO) flzz(O) 



Ml 



si n^^iMl) s„ 2nzy{Ml) 



Ml 



Ml 



(A3a) 
(A3b) 

(A3c) 



The hats denote self-energies renormalized using the on- 
shell renormalization condition: 



^ Rc YlzziMl) 

= nz7(o) = n^^(o) = 0. (A4) 



Explicit formulae for the renormalized self-energies in 
terms of the unrenormahzed ones are 16911 : 



^2 



Mi 



3-iU-n' (0) 



where 



R 



Uzz{Ml) Uww{M^) ^s.^Uz^{0) 



Ml 



Ml 



(A5a) 

(A5b) 

(A5c) 
(A5d) 

(A6) 



Substituting (|A5|) into (jA3[) and taking into account that 
nz7(0) = we find: 



4s^ 



nzz(O) - TlzziMl) + li^yiMl) 



'-^z^iMl) 



T 



TlzziO) _ UwwiO) 

nww{o)-n 



Ml 



(A7a) 
(A7b) 



U = 4s^ 



WW 



nzz(o) 



nzz(M|) , slU^,{Ml) 



Ml 



Ml 



(A7c) 



Appendix B: Calculation of S,T,U 

To evaluate the contribution of n Majorana neutrinos 
to the S, T, U parameters defined in Appendix|X]we need 
to calculate the self-energies entering (|A7|) . Only charged 
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FIG. 15: Contribution of the charged leptons to 11^^ and IIz^ 
at one-loop level. 

leptons contribute to 1X27 and II^^, see Fig.[Tni The 
resulting self-energies are the same as in the SM: 



Tl^iq') = 



mlBoiq^ml,ml)] , (Bla) 
(Bib) 



1) 



where m„ denote the masses of the charged leptons. To 
shorten the notation in (IBll) we have introduced 



+ q^[Bi{q\ml,ml) + B2i{q^,nil,ml)\ , (B2) 

where Bn, B21 and B22 arc the usual one-loop func- 
tions [lOl, = 4 — 2e and e — > 0. The one-loop con- 
tribution of the charged leptons and neutrinos to IIvi/w 
is presented in Fig.[Tni Note that due to the nonzero 




FIG. 16: Contribution of the charged leptons and neutrinos 

to Ilvi^ty. 

active-sterile mixing the heavy neutrinos can also run in 
the loop. In position space the resulting self-energy is 
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given by: 



2sl ^ 



lUo 



^io.{x)^^PLy\x)v,{y)^''PLt''{y)\). (B3) 



Since we deal with Majorana fcrmions, to evaluate (jB3[) 
it is convenient to use the two component notation for the 
spinors, see Appendix[C]for more details. The self-energy 
then takes the form: 



X {T{xAx)ani{x)ll{y)-o^x\y) 

+ ll{x)-a^^{x)x^{y)-o''ll{v)\). (B4) 

where Ij^ denotes left-handed component of the charged- 
lepton field in the two-component notation. To evaluate 
(|B4p we need to consider all possible contractions of the 
field operators. Using the Fourier-representation of the 
propagators we obtain: 



24 



tr \$-i^fY\ 
(jp' — -I- ie) (fc^ — + ie) 



(B5) 



Note that since only the left-handed component of the 
charged field runs in the loop, only 'kinetic' contraction 
is possible. This is reflected by the fact that the nu- 
merator of (|B5p contains only the momenta ^ and f of 
the intermediate states. Taking the trace and using the 
definitions of the one- loop functions we finally arrive at: 



- P{^.rn\, TO^)] , (B6) 



where Q has been defined above (|B2[) and 

P{q^ml ml) = 2B2i{q\ml, ml) + 2Bi(g2, m?, m^) . 

For the Z boson self-energy we need to consider two di- 
agrams with the neutral or charged states propagating 
in the loop, see Fig.[T71 In this case both 'kinetic' and 
'mass' contractions are possible. In complete analogy to 
()B6p the kinetic contraction of the two neutrino lines re- 
sults in: 



nrz(i)('z^) 



^ ^ tj 



^[9ti.vQ{q^ ,ml,m)) - q^qy P{q^ ,mj,m])] . (B7) 





FIG. 17: Contribution of the charged leptons and neutrinos 
to lizz- 



Since the intermediate neutrinos are Majorana particles 
the 'mass' contraction is also possible: 



^^^^EEuLu.,ut,u,, 



X [gf^ymimjBo{q^,mf,mj 



(B8) 



This contribution vanishes for vanishing Majorana mass 
as it should. Note also that the flavor structures of (jBSp 
are slightly different compared to (jB7|) . 

Since the left- and right-handed charged leptons couple 
to Z with the strengths 1 — 2s^, and 2s^ respectively, the 
contribution of the kinetic contraction reads: 



;[(i-24f + (24)^1 



^[g^i^Qiq^ ,ml,ml) ~ q^q^ P{q^ ,ml,ml)] . (B9) 



Since charged leptons are Dirac particles the mass term 
appears only for contractions of the left- and right- 
handed components and is therefore proportional to a 
product of the two couplings: 



(a'^) - — 



^44(1-24) 
x[g^,,mlBo{q^,ml,ml)]. (BIO) 



The total contribution to the Z boson self-energy is given 
by the sum of Eqs. (|B7| - ((BT0| . 

Substituting the terms proportional to g^'^ of 
and 11^2 into the deflnitions of the S, T, U parameters 
we obtain Eqs. ([TSj) . (fT8|) and ([T9| . Since the one- loop 
integrals are divergent, 

Q^'^q^ml ml) = e-i(m?/2 + ml/2 - q^i) , (Blla) 



B^^^q',mi,mi) 



(Bllb) 



each of the terms in these expressions is divergent as well. 
However, their combinations in Eqs. (fTSj) . (fT8| and (fT9|) 
arc flnitc. This can be checked explicitly by using the uni- 
tarity of the full mixing matrix U as well as the relation 



{mL)ap = Z]j=i" Uaim^U^ = 0, where a,;9 G {e,/i,T}. 
The latter reflects the fact that for a type-I see-saw the 
upper-left corner of the mass matrix is zero in the fla- 
vor basis. Since the divergences and the scale /i in the 
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finite parts of the loop integrals always appear in the 
same combination, + In/i^, the cancellation of the 
divergences implies that the /i-dependence drops out as 
well. In other words, the S, T, U parameters depend only 
on physical quantities like the couplings, masses and mo- 
mentum transfer. This is a consequence of the fact that 
they are defined in terms of the on-shell renormalized 
self- energies, see Eq. (jA3p . 



Appendix C: Two-component notation 

To evaluate the contribution of the Majorana fcrmions 
to the self-energies it is convenient to use the two- 
component spinor notation. In terms of the two- 
component spinors the four-component Dirac spinor and 
its Dirac-conjugate read [TSj : 




(CI) 



with £, = cl and £, = cr for the charged leptons. For 
Majorana fcrmions the two spinors in (jCl[) are not inde- 
pendent, ^ = Xi s-iid therefore the number of degrees of 
freedom is reduced from four to two: 



*M = 



X'' 



{x'.Xp): 



(C2) 



with X = ^ for the Majorana neutrino. The Dirac matri- 
ces can be written in the form 



(C3) 



where ct" = ct" = 1 and ct^ = — cr^ with cr^ being the 
Pauli matrices for /i = 1,2,3. Using this representation 
and the known formula for traces of the Dirac matrices 
we find: 





tr[7^7'^] = tr^CT''] + tY[a^'(J''] , 



(C4a) 



tr[7^7'^7''7^] = tTla^'a^aPa^] + tila^" a" a" a^] . (C4b) 



Using (jC2[) and the known form of the Feynman propa- 
gator of a Dirac field we can now infer the form of the 
propagator of a Majorana field: 

(27r)4 p2 _ „,2 -^iA p^^i^Jj m 



(C5) 



The diagonal components of (jCSP describe contractions 
of the field with itself and reflect the Majorana nature of 
the fleld. 



To use the above formulae we need to rewrite La- 
grangian ^ in terms of the two-component spinors: 



/2s„ 



where a are the flavor and /3, /3 the spinor indices. As can 
be inferred from (|C6[) . the contribution to the Z boson 
self-energy is proportional to 

{T[D^{x)a'^^f^^v,{x)Df,{y)a''^^^v,{y)]), (C7) 

where we have suppressed the generation indices to 
shorten the notation. Now we need to use Wick's the- 
orem and find all possible contractions. From (|C5[) it 
follows that for Majorana fermions there are two possi- 
bilities. The first is the contraction of the field with its 
conjugate. This gives rise to 

cx -tr[p • aa^k ■ era"] ^ -itr[^7^|^7''] . (C8) 

The second contraction possible only for Majorana 
fcrmions is the contraction of the field with itself. It 
gives rise to 

- {TD^ {x)D^ (y)) . (Tz., (y),,^ {x))a'^'^^ 

oc -mimjtr[cr^g-''] = - im,;mj tr[7''7''] , (C9) 

where we have used the identity 



X/3<^ X7 = -x'^;^ X^ 



(CIO) 



Collecting the two contributions and taking the traces of 
the Dirac matrices we find that (jC7P is proportional to: 



-2 {p^'k'' + p^k - (pk) gf"" + rmm^g'"') . (Cll) 

The contribution to the self-energy of the W boson is 
proportional to 

(r[i>^(x)a''^'^^ei^(x)e-i^(y)a''''^V,(y)]) . (C12) 

Because of the Dirac nature of the charged leptons only 
one, namely the 'kinetic' contraction is possible in this 
case. It results in an expression identical to 



Appendix D: Metropolis Algorithm 

In this work we minimize the function in a multidi- 
mensional parameter space. Due to the existence of strict 
constraints on the observables, the gradient minimization 
methods are not efficient for finding global minima. In- 
stead a statistical method is applied. In the first step we 
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perform a random scan in the nine dimensional param- 
eter space. The scan parameters are the three complex 
angles of the arbitrary orthogonal matrix in Eq. (|20|) and 
the three masses of the sterile neutrinos. The masses are 
chosen on a logarithmic scale to cover most efficiently a 
broad mass range. 

If a point of the parameter space chosen in this way 
satisfies the — >■ 67 and QvjSp constraints, a second step 
is performed. In the second step we calculate the full 
and perform a local minimization. For this purpose we 
utilize the Metropolis algorithm which is usually taken 
to simulate phase transitions in the Ising model. First a 
fictive parameter (which we call temperature T) is intro- 
duced. Then parameters of the potential good-fit point 
are changed randomly by a small amount. The Xncw 
computed for the new point and the Boltzmann function 
B{T,Ax^) with the difference is evaluated: 

B(r,Ax^) =exp f-1^^2id_^2^y (Dl) 



Furthermore, following the equal probability distribu- 
tion, a random variable < x < 1 is generated. If 
Xncw < Xoid or B{T,Ax'^) > X the new point is cho- 
sen as the new starting point. If B{T, Ax^) < x the 
new point is discarded. This process is repeated at a 
temperature 71 until a quasi-cquilibrium is reached and 
no large changes in occur. Then the temperature is 
decreased to 72 < 7i and the process is repeated until 
the quasi-equilibrium is reached at the new temperature. 
This iterations continue until the effective temperature 
is zero. 

This method proves to be more efficient than a gradi- 
ent method for finding global minima in the considered 
case since the system is highly constrained and thus the 
parameter space has a very non-trivial topology. Using 
the finite effective temperature approach we are able to 
scan a larger part of the parameter space before the sys- 
tem settles in a minimum when the temperature drops 
to zero. Therefore, with a higher probability, the found 
minimum is the global one. 
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